In [1]:
# Standard setup block for running Python code
from __future__ import division, print_function
import os
if os.path.split(os.getcwd())[-1] == "Lab notebooks":
os.chdir("../../")
print("Moved to experiment root directory")
from Modules.processing import *
from Modules.plotting import *
plt.style.use("Config/plotstyle.mplstyle")
%matplotlib inline
Got to the lab about 3:30 PM.
Homed turbine, tow, y-, and z-axes.
Tank depths is 2.44 m -- good.
Plenty of seeding particles have risen to the surface.
3:35 PM -- Did two dummy tows at 1 m/s to mix up seeding in the tank.
Used squeegee to mix seeding particle from surface. Vectrino SNR looks very good--up around 14 dB.
3:47 PM -- Re-homed y- and z-axes to be safe. Everything looking good.
3:54 PM -- Restarting Wake-1.0-0.5.
4:20 PM -- Serial device server is acting weird. Reset and things are okay again.
4:32 PM -- Continuing Wake-1.0-0.5 from run 9.
4:55 PM -- Noticed that halfway through run 11 Vectrino stopped recording. Need to redo that run. Redid the run, but the data was actually fine. Was being over-filtered due to some high turbulence intensity from tip vortex shedding.
8:04 PM -- Completed Wake-1.0-0.5. z-axis RPOS at lower limit it -0.03600 m. y is -1.52549 m.
8:11 PM -- Starting Wake-1.0-0.625.
9:40 PM -- Finished Wake-1.0-0.625 up to run 16. Ending for today due to incoming blizzard. y-axis at negative limit RPOS = -1.52574 m.
In [14]:
Section("Wake-1.0-0.625").process(nproc=1, nruns="new")
In [15]:
w = WakeProfile(1.0, 0.625)
w.plot("mean_u", "-o", preliminary=True)
In [4]:
r = Run("Wake-1.0-0.5", 15)
plt.plot(r.time_vec, r.u)
plt.show()
In [5]:
length = 0.5
diameter = 2.9*0.0254 # inches to m
area = length*diameter
cd = 1.2
r_force = 0.25
rho = 1000.0
u = 1.0
tsr = 5.5
f = 0.5*rho*area*cd*(tsr/2*u)**2
torque = f*r_force
print(f, torque)
In [11]:
# Setting things up
setup_time_sec = 3600
n_setup_times = 4
# Seconds per tow at each speed
seconds = {0.8 : 90 + 200,
1.0 : 84 + 240,
1.2 : 80 + 300}
# Wake runs
runs_per_wake = 43
nwake = 1 + (43 - 14)/43
wake_speed = 1.0
# Runs for strut torque and strut covers
strut_covers_runs = 22
strut_torque_runs = strut_covers_runs*2
# Tare drag and torque
tare_drag_runs = 14
tare_torque_runs = 30
tare_drag_dur = 200
tare_torque_dur = 120
total_time_sec = n_setup_times*setup_time_sec + runs_per_wake*nwake*seconds[wake_speed] \
+ strut_covers_runs*seconds[wake_speed] + strut_torque_runs*tare_torque_dur \
+ tare_drag_runs*tare_drag_dur + tare_torque_runs*tare_torque_dur
total_hours = total_time_sec/3600
hours_per_day = 10
total_days = total_hours/hours_per_day
print(total_hours, total_days)